Theory of Dephasing by External Perturbation in Open Quantum Dots 
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We propose a random matrix theory describing the influence of a time dependent external field 
on the average magnetoresistance of open quantum dots. The effect is taken into account in all 
orders of perturbation theory, and the result is applicable to both weak and strong external fields. 
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It is well established that the anomalous magnetore- 
sistance of bulk disordered systems is governed by the 
weak localization ( WL) . Being an interference phe- 
nomenon, the WL is extremely sensitive to inelastic pro- 
cesses which are commonly referred to as dephasing. 

Recently, another object for the studying the quan- 
tum effects appeared - ballistic quantum dots |^ . In the 
absence of inelastic processes, the transport properties 
of the dots are well described within the Random Ma- 
trix Theory (RMT) ||^. The magnetoresistance within 
this theory manifests itself as a crossover between two 
universal ensembles (orthogonal and unitary), and the 
strength of the magnetic field defines a position on that 
crossover. This approach per se does not include de- 
phasing, and the dephasing processes were considered on 
a phenomenological basis The relation of this phe- 
nomenological description used to fit the data of Ref. |^] 
to microscopic mechanisms of dephasing is still an open 
question. 

In this Letter, we propose a Random Matrix - like the- 
ory of the magnetoresistance affected by an external time 
dependent perturbation. We will be able to find both the 
amplitude and frequency dependence of the magnetore- 
sistance using only one unknown parameter. This param- 
eter can be related to the correlator of the level velocities 
due to the same perturbation at zero frequency and, thus, 
in principle can be determined by an independent exper- 
iment. After the strength of the potential is normalized 
by this parameter, all the results become universal. From 
the experimental point of view, we have in mind chang- 
ing shape of a quantum dot by applying an external ac- 
bias. 

Before we proceed, let us mention that the effect stud- 
ied in the present paper is similar to that of Ref. 
where it was shown that uniform ac- electric filed sup- 
presses the weak localization correction to the conduc- 
tivity of a disordered wire (experimentally it was studied 
in Ref. Q). However the results of are not directly ap- 
plicable to the quantum dots with size L so small that the 
Thouless energy ~ h/r^j-g is much greater than other 
energy scales (such as the dephasing or escape rates) of 
the problem (here r^rg is the characteristic time for the 
classical particle to cover all of the available phase space) . 



On the other hand in this limit one can use the RMT 
to study the conductance of the system, see [|l0| . All cor- 
rections to the RMT are small as Nch/gdou 9dot = Et/5i 
and 5i is the mean level spacing. We consider the WL 
correction to the conductance of quantum dots with the 
large number N^h of open channels. In this approxima- 
tion we neglect the effect of interaction on the conduc- 
tance which are small as XjNl-f-^ while the weak lo- 
calization correction to the conductance is proportional 
to 1/Nch- The same condition also allows us to use con- 



ventional diagrammatic technique 1 12| to perform the en- 
semble average. 

In general, the Hamiltonian of the system can be rep- 
resented as, see llToll: 



H = Hr 



Hr+H 



LD, 



(1) 



where Hu is the Hamiltonian of electrons in the dot, 
which is determined by AI x AI matrix 



M 



(2) 



where the thermodynamic limit M ^ oo is assumed. We 
consider the case, when -ff„„i is a time dependent random 
matrix in the form: 



(3) 



Here the time independent part of the Hamiltonian Tinm 
is a random realization of M x M matrix, which obeys 
the correlation function 



(4) 



where A = MiSi/n)'^ and A' = A(l - gh/4:M), gu defines 
crossover from the orthogonal {gh — 0) to the unitary 
{gu — 4M) ensemble. Parameter gu has a meaning of the 
dephasing rate by the external magnetic field in units of 
the level spacing 5i, see |§,|lO|. It can be estimated as 
gh — gdot i^/^o)^ where $ is the magnetic flux through 
the dot and $o = hc/e is the flux quantum. The time de- 
pendent perturbation is described by symmetric M x M 
matrices Vnm and a function of time (p(t). 

The coupling between the dot and the leads is 



1 



H 



LD 



(5) 



where V'n correspond to the states of the dot, ipa{k) de- 
notes different electron states in the leads, and momen- 
tum k labels continuous state in each channel a. For the 
dot connected with two leads by Ni and Nr channels re- 
spectively, we denote the left lead channels by 1 < a < A^; 
and the right channels hy Ni + 1 < a < Nch, where 
Nch — Ni + Nr- The electron spectrum in the leads near 
Fermi surface can be linearized: 



HL^vpj2k^i{k)^^{k), 



(6) 



where vp — l/inv is the Fermi velocity and v is the 
density of states at the Fermi surface. 

The coupling constants Wna in Eq.(^)are defined as 



W - / if = a < Nch, 

- Y ^2 J, I'o, otherwise, ^ > 

where ta determines the dimensionless conductance of 
each lead (in units of 2e'^/h) according to 



9r 



^"'^ Ai i* 



(8) 



and \ta\ < 1. The factor in Eq. (Q) is chosen so that the 
ensemble average scattering matrix S^p of a dot with 
fully open channels (<q — 1) is zero. More complicated 
structure of W can be always reduced to the form by 
suitable rotations. 

For the system described above the scattering matrix 
S has the form: 

Sa.0{t,t') = 1 - 2mi^Wl^Gnm(t,t')W^(3, (9) 

and the Green's function Gnm{t,t') is the solution to: 
0^ - H{t) + iTTi^WW^^ G{t, t') = 5{t - t'), (10) 

where matrices H and W are comprised by their elements 
(^ and (0) respectively. 

The averaged dimensionless dc - conductance of the 
dot is determined in terms of the scattering matrix of the 
system in the linear response theory by (see, e.g. 0): 



9 = 



(11) 



where (. . .) stands for both ensemble and time averages. 
We also introduced notation for the projector on the left 
lead, f/, which is a diagonal Nch x Nch matrix with the 
first Ni diagonal elements equal to unity, and the other 
diagonal elements equal to zero, and fv- = I — fi. 



We perform calculations of the average conductance 
keeping the leading terms in 1/M. The diagrammatic 
technique is somewhat similar to that developed for bulk 
metals [T^ ], where the small parameter is l/epTimp with 
ep being the Fermi energy and Ti^p being the elastic 
mean free time. 
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FIG. 1. (a) Diagrams for the ensemble average Green's 
function. The second term in the self-energy includes an in- 
tersection of dashed lines and it is small as 1/M. (b) The rep- 
resentation of the conductance in the form of Eq. (|l3) forbids 
the renormalization of vertices J from Eq. by disorder. 

First let us find the ensemble average Green's function 
(G(^)). One can see that (G^-^)) is diagonal, (G^^(e)) = 
SnmGi[^\e) Using the self consistency equation for the 
Green's function, Fig. 1 (a), we find 



1 



1 + f t* 



-I- ie 



AM 



n < Nch, 

n > Nch- 
(12) 



Here we introduced the dimensionless energy e measured 
in units of a/A/4M — Si/2tt. We expand these Green's 
functions in e/M and {gi + gr) / M , since only those terms 
survive the thermodynamic limit M —> oo. For the same 
reason, the expression for G^^' with n < Nch neglects 
such terms at all because the contribution of those ele- 
ments to the final answer is already small as Nch/M. 

To simplify further manipulations, we rearrange 
Eq. (|ll| ) in the following form 



g{t)^(^ J dt'tY 

— OO 

^ Nrgf + Nigl 

{gi + 9rY 

91,: 



JlS{t,t')JrSHt,t'] 



Jl.r — "^Lr 



9l+9r 



■I. 



(13) 



(14) 



Equation (|1^) immediately follows from Eq. ( |llD and the 
unitarity of the S- matrix 55^ = 1. The calculations of 
the conductance in the form of Eq. (^) are significantly 
simpler since the vertices (^4|) are not dressed by dashed 
lines, see Fig. 1(b). This trick is similar to the calculation 
of the conductivity of disordered bulk systems in terms of 



2 



the current-current instead of density-density correlation 
function, see Refs. 

Now we substitute the scattering matrix defined by 
Eq. d) to Eq. (|l|). To the leading order in l/{gi+ gr) 
one can average S -matrices independently with the help 
of Eq. (|l2) and obtain the classical conductance 



9l9r 



m 



(15) 



In particular, for the dot with fully open channels {ta — 
1), the averaged S- matrix vanishes {{S) = 0) and the last 
term of Eq. ( p^ ) gives the known result gd ~ NiNr/Nch, 
since in this case gi^r — N^r- 
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FIG. 2. (a) The diagram for the WL correction to the 
conductance, (b) The diagram equation for the Cooperon. 

The first order correction in l/{gi + gr) to Eq. ( |l5| ) 
is given by the diagram in Fig. 2(a). It represents the 
WL correction to the conductance, and has the analytic 
expression 



figl + /r.9f [ CjdT 



[gi + 9rY J 27r 





2dTC{T,T,-T), (16) 



where formfactors fi^r are given by 



16(iaC)' 



fr 



y - 

^ ^ (1 



(l + M*)^ 



The Cooperon C is defined by Fig. 2(b): 

'|- + /C(T,r))c(T,T,T')=5(r-r'), 



(17) 



(18) 



where time is measured in units of inverse level spacing 
2ti/5i and the "Hamiltonian" for the Cooperon is 

/C(r, r) = 5, + TT^Co b(T + t/2) - ^{T - T/2)f , (19) 



with g^ characterizes the total dephasing due to the es- 
cape and the magnetic field, = 5z + ffr + and we 
chose if{t) = cosujt to describe the time dependence of 
the perturbation. 

The only unknown parameter, Co, in Eq. ( p^ ) depends 
on the strength of the perturbation. In terms of the orig- 
inal Hamiltonian (^, it is defined as 



(20) 



where we used the fact that the matrix V is symmet- 
ric. This parameter is also related to the typical value 
of the level velocities, which characterizes the evolution 
of energy levels e^{X) under the action of the external 
perturbation XV, see |l3[: 



dX 



\dx/ 



(21) 



Since all other responses (e.g. parametric dependence 
of the conductance of the dot) are expressed in terms 
of universal functions of the same parameter Co [ p^ , 
it can be found from independent measurements. For 
not very realistic case of homogeneous electric field E in- 
troduced into the dot of linear size L, one can estimate 
Co — (eEL)^ /{EtSi). It is important to emphasize that 
the homogeneous shift of all levels does not affect the 
magnetoresistance and that is why the average level ve- 
locity (de^/dX) is not relevant. 

In the absence of the time dependent perturbation 
= 0, one obtains from Eqs. (|l|)-(|l9|): 



figr + frgf 



^^""^ {gi+gr^g* 



(22) 



The solution to Eq. ( |18| ) gives the weak localization cor- 
rection to the conductance Ag^i in the presence of the 
time dependent field. It can be expressed in terms of the 
unperturbed correction (p^) as 



Agwi 



g*5i 



g* 



(23) 



where dimensionless function F{y, z) is given by 



Jo 



smxy 
= x -. 24) 

y 



Here Io(0 is the modified Bessel function. Some curves 
for this function are plotted in Fig. 3. 

Equations ( |2^ ) - (24) are the main results of our pa- 
per. They give the universal description of the effect 
of the external field on the weak localization correction. 
Below we will discuss different asymptotic regimes and 
compare them with the results for bulk systems ^ . 

For weak external field z <C max(l, y"^) we find 



3 



As 



wl 



= 1 - 



5W 



(25) 



In this regime the correction is quadratic in the frequency 
for slowly oscillating field, similarly to the bulk system re- 
sult at w smaller than the dephasing rate I/t^. However, 
the frequency dependence saturates at large frequency. It 
is different from the result for bulk systems, where a char- 
acteristic spatial scale shrinks as whereas in our 
case it is determined by the size of the dot. 
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FIG. 3. Representative curves of F(y, z) as a function of 
z for two values of y. It decreases linearly with z at small val- 
ues of z. The inset shows the j/— dependence of the function 
F{y, z) for two values of z. It decreases quadratically in y at 
small values of y and saturates at larger y. 

In the opposite limit of strong external field z 3> 
max(l,?;~^) we have to consider separately the cases of 
fast, y ^ I, and slow, y <^ 1 field oscillations. In the 
first case we find 



^9wi 
A.i? 



9* 



(26) 



The linear dependence of the quantum correction on 
l/x/Co is similar to that for the bulk system. Contrary 
to the bulk systems, the result does not depend on the 
frequency for reasons we have already discussed. 

In the case of slow field y ^ 1, but still zy"^ ^ 1 (strong 
field) we obtain 



Ag^; _ r(l/6) ( 25lgl 
A^^ 



7rr(5/6) VSCow^ 



1/3 



(27) 



i.e., the dependences both on the amplitude and fre- 
quency are different from the bulk case. 

In conclusion, we proposed a random matrix theory 
describing influence of time dependent external field on 
the average magnetoresistance of open quantum dots. 



This dependence can be recast in the form of the uni- 
versal function Eq.(p^ of one fitting parameter Eq.([2l|) 
which can be fixed by an independent experiment. The 
results can not be described by a simple replacement 
g* ^ g* + 14,. Finally, we mention that thermal fluctua- 
tions of the gate potentials may induce the dephasing by 
virtue of the mechanism considered here. However, the 
spectral density of such fluctuations is model dependent 
which makes quantitative predictions hardly possible. 
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